Figure Captions 



Figure 1: The thick- walled bubble profile. 
Figure 2: The spectrum of created particles. 
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■ Abstract 

O ! 

Q . The spectrum of created particles during the tunnehng process, leading to the decay of a 

CN \ false vacuum state, is studied numerically in the thick-wall approximation. It is shown that in 

K*- ' this case the particle production is very intensive for small momenta. The number of created 

^ ^ , particles is nearly constant n{p) ~ 1 for 4 < p < 500. 

\ 

in ' 

^ 1 Introduction 

I> ' In the semiclassical approximation, false vacuum decay through the barrier penetration is described 
^ by the nontrivial 0(4) symmetric solution of Euclidean (imaginary-time) equation of motion with 
psj ! the boundary condition that the field asymptotically approaches its false vacuum value, This 
^ \ solution is referred to as the bounce one. The bounce looks like a four-dimensional spherical bubble 
' of true vacuum embedded in a false vacuum background. The bubble evolution after nucleation is 
O obtained by the analytic continuation of the bounce solution to the Minkowskian time, r — it. There 
• exist two limiting cases, referred to as the thin and thick wall approximations respectively, when the 
^ \ theory of bubble formation is greatly simplified, |2] . Thin wall approximation is realized when 
Q \ the energy- density difference between the true and false vacuum is much smaller than the height of 
'-^ ■ the barrier. The particle production in this case has been studied in jHj. The analysis presented in 
that paper shows that the particle production is strongly suppressed in the thin-wall approximation. 
^ It is of interest to explicate whether this is the general feature for the created particles spectrum 
during the false vacuum tunneling when the background space-time curvature is neglected (without 
gravity). Throughout this paper we shall consider another limiting case, when the energy-density 
difference between the true and false vacuum is much grater than the height of the barrier. The 
paper is organized as follows. In Sec. 2 we briefly describe the technique for calculation of the particle 
spectrum. In Sec. 3, based on the calculational procedure described in Sec. 2, we evaluate the particle 
spectrum for the thick-walled bubble. Sec. 4 is devoted to the summary and concluding remarks. 



J3 



2 Technique of calculation 

For the calculation of particle spectrum we use the standard technique developed in jHEj. Through- 
out this paper the metric signature in Minkowski space-time is (— , +, +, +) and units U = c = 1 are 
used. We consider a scalar fleld deflned by the Lagrange density 



with the potential 



C = -\{d,ct>f-U{ct>), (1) 



U{<P) = ^02(0 _ 2af - e(0V2a3 - 30Vl6a^), (2) 



where A, a and e are positive parameters. The potential (j2)) has a local minimum at (j) 
0, a global minimum at (pt = 2a, U{2a) = — e, and a local maximum at (pt^p = 
3e), U{(f)top) = l28X^a^\2Xa^ + e)/{8Xa'^ + 3ef. The 0(4) symmetric bounce solution 
Euclidean equation of motion 

^ + = U'((j)) 
dp^ p dp 



= 0, U{0) = 
UaV(8Aa^ + 
<h satisfies the 

(3) 



with the boundary conditions 0(oo) = 0, dip /dp = at p = 0, where p = a/x ^ + and the 
prime denotes differentiation with respect to 0. Since the bounce solution is 0(3,1) invariant in 
the Minkowskian space-time and 0(4) invariant in the Euclidean region it is convenient to use the 
following coordinate systems. The coordinate system used in the Euclidean region is (p, x, 6', ip) where 
{6, if) are usual angle coordinates on two-dimensional sphere and (p, x) are related to r = \x\ and r 
as follows 



r = psin(x), r = -pcos(x), 
< X < 7r/2, < p < oo. 

The coordinate system in the Minkowskian region may by obtained by replacement (p, x) 
{-iphi, -iXM), which yields 

r = pm sinh(xAf), t = pm cosh(xM), 
< XAf < oo, < Pm < oo. 

The equation governing the fluctuation fleld $ reads 



(4) 



M 



Pm Pm 



$ = 0, 



(5) 



(6) 



where denotes the partial derivative with respect to pm and L is the Laplacian operator on three- 
dimensional unit hyperboloid. Expanding $ in terms of harmonic functions on the three-dimensional 
hyperboloid Ypim, —L^Ypim = (1 +p^)l^«m, the Eq.® for the mode function takes the form 

{p' + 1/4) 



(9^ -u 



pIi 



+ U"{ 



^ = 0. 



Making use of analytic continuation of ip to the Euclidean region by replacement pm 
the equation for the mode function in the under-barrier region 

(p2 + 1/4) 



p- 



-U"{ 



0, 



(7) 

ip one obtains 
(8) 



where ipE is connected with ip by the asymptotic boundary condition iPe{.—Wm) = "^{pm) at pm 0. 
At r = — oo the fleld is in false vacuum state and correspondingly the fluctuation fleld satisfles 
the vanishing boundary condition tpE ^ when p ^ oo. Since U"{(j)b) —>■ 4a^A when p oo the 
asymptotic solution of Eq. ((7j) will have the form 



M 



Ci{p)y/p^ Hip\2ay/XpM) + C2{p)^/pli H\f{2a\/\pM) 



(9) 



where Hj;p\x) and H-^\x) are the Hankel functions of the flrst and second kinds, respectively. The 



spectrum of created particles n{p) is given by jHE] 

n{p) 



-2pTT 



,|ci(p)/c2(p)p - e-2p^| 
For detailed description of the general formalism see |31E1- 



(10) 



3 Particle spectrum 



For the thick-wall case we specify the potential parameters as follows a = 1 , X = 0.5 , e = 100. For 
that values of parameters the ratio U{(j)top)/^ = 0.5752022890e — 6. The corresponding bounce 
solution obtained by the program code odex with local tolerance l.OD — 17 is shown in Fig.l. 
The detailed description of this program code is given in 6J. Taking into account that 0fe(lO) = 
0.3453646895630613e — 08 and using the program package Maple 7 one finds that the solution of 
Eq.(jH)) for p > 10 satisfying the vanishing boundary condition when p cxd is given by 



i^E = a3y/pKip{y/2p), 



where Kip{x 



is the modified Bessel function of the second kind. In order to make the analytic 
continuation at p = we construct the analytic solution to Eq.(jH)) in the vicinity of p = 0. For this 
purpose, the power series expansion of (pi, about p = is used for evaluating of U"{(j)b) 



-29.34386061 + 50.39180321p2 - 52.71178464p^ - 21.76409860p^ + 0(p^ 



(12) 



Since we want to use the expansion ()12j) for evaluating the solution of Eq.(jH)) in the region < p < 
0.001 we omit the terms higher than quadratic. Inserting this expression into Eq.(jHJ and using the 
program package Maple 7 one obtains the following solution 



ai 



a2 



M(1.033421066, 0.5 zp, 7.09871842p^) + —1^(1.033421066, O.bip, 7.09871842p^), (13) 



Vp VP 

where M, W denote the Whittaker functions. For pM > 10 the solution of Eq.((Z|) is given by 

^(Pm) = Clip). 



{V2pm) + C2{p)y/PM H^p{V2 Pm)- 



(14) 



Since the particle number depends on the ratio \ci{p)/c2{p)\ the multiplier in Eq. (fTT| may be 
fixed arbitrarily. The Eq. ljllj) and its derivative give the initial conditions to Eq.® at p = 10. Using 
these initial values the solution of Eq.(jHJ in the region 0.001 < p < 10 is constructed numerically by 
the code odex with local tolerance l.OZ^ — 17. Matching this solution with Eq. ()13|) at p = 0.001 one 
finds the coefficients ai, a2- Making use of analytic continuation p —>■ —ip the Eq. p3|) now gives the 
initial conditions to Eq.Q at pm = 0.001. With these initial conditions the Eq.Q is solved again 
numerically up to pm = 10 by the odex with local tolerance l.OD — 17. Matching this solution with 
Eq.dlll) and using Eq.()10p one evaluates the particle number for given p. For evaluating of special 
functions and their derivatives at the matching points we have used the program package Maple 7. 
The spectrum obtained in this way is shown in Fig.2. As p increases from to 4 the number of 
created particles decreases from 23714.66515 to 1.000916176 and becomes nearly constant n{p) ^ 1 
for 4 < p < 500. For the sake of clarity the behavior of n{p) for small values of momentum is shown 
in the following table. 



p 


0.25 


0.5 


0.75 


1 


1.25 


n 


15524.29139 


6622.761821 


2612.820881 


1008.019999 


378.9347624 


P 


1.5 


1.75 


2 


2.25 


2.5 


n 


138.1763281 


49.01973850 


17.20855685 


6.281386004 


2.664255183 



4 Summary 

Using the standard technique |31IS] we have shown that, in contrast to the thin-wall case, the particle 
production in the thick-wall approximation is very intensive. The number of created particles is 



especially large at small values of momentum, decreases sufficiently fast as p increases from to 4 
and becomes nearly constant n{p) ^ 1 for 4 < p < 500. To date we have been unable to evaluate the 
behavior of n{p) for asymptotically large values of p. In the thin- wall case the particle production is 
strongly suppressed in general, it is nearly constant for small momenta and behaves as exp(— 2p7r) 
for large values of momentum j3] . As we observe the particle production has a quite different nature 
in the thick-wall case. A few remarks are in order. This result is not strictly correct since we do 
not take into account the constraint on the fluctuation field that arises due to fact that the proper 
fluctuation field associated with the tunneling is the transverse part of total fluctuation field with 
respect to the bounce solution [7 . In the application of this constraint is illustrated for the 
spatially homogenous tunneling. It was shown that, due to this constraint, the spatially homogenous 
tunneling does not allow the particle production with zero momentum. If such a restriction will 
not reduce the particle number significantly in the thick-wall case, one has to consider the back 
reaction of particle production on the tunneling. On the other hand the analysis of a concrete 
thick- wall model, when the quartic term is ignored in Eq.(0), shows that the activation rate at zero 
temperature makes a dominant contribution to the vacuum decay in comparison with the tunneling 
rate jHj. Unfortunately we do not know whether this is valid in the thick- wall approximation in 
general and if this is the case how does the zero temperature activation rate depend on a^A/e. 
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